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ABSTRACT 

We report the serendipitous detection by the XMM-Newton X-ray Observa- 
tory of an X-ray source at the position of the Type I (He- and N-rich) bipolar 
planetary nebula Hb 5. The Hb 5 X-ray source appears marginally resolved. 
While the small number of total counts (~ 170) and significant off-axis angle of 
the X-ray source (~ 7.8') precludes a definitive spatial analysis, the morphol- 
ogy of the X-ray emission appears to trace the brightest features seen in optical 
images of Hb 5. The X-ray spectrum is indicative of a thermal plasma at a tem- 
perature between 2.4 and 3.7 MK and appears to display strong Neon emission. 
The inferred X-ray luminosity is Lx = 1.5 x 10 32 ergs s _1 . These results suggest 
that the detected X-ray emission is dominated by shock-heated gas in the bipolar 
nebula, although we cannot rule out the presence of a point-like component at 
the position of the central star. The implications for and correspondence with 
current models of shock-heated gas in planetary nebulae is discussed. 
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Subject headings: planetary nebulae: individual (PN Hb 5) 
ual(PN Hb 5) 



X-rays: individ- 



1. Introduction 

The shaping and evolution of planetary nebulae (PNe) has been an active field of study 
for a few decades an d is continually spurre d on by new observations that do not fit the existing 
theoretical models (IBalick fc Frankll2002l . and references therein). Of particular note are the 
bipolar PNe, with their large lobes and narrow waists, many showing residual evidence of a 



spher ical AGB wind in the form of nested rings centered on the PN central star (ICorradi et al. 
20041 ). The implied transition from seemingly spherical mass loss to profoundly bipolar 



outflow remains an unsolved mystery. Possible shaping and evo lutionary mechan isms have 
been developed ranging from the interacting stellar winds model (IKwok et al.lll978l). wherein 



a hot fast stellar w ind shocks a slower moving AGB wind, to collimated flows (jSokerl 12004 



Akashi et al.ll2008l) possibly launched a nd shaped by a magnetocentrifugal force fueled by 
an AGB dynamo (IBlackman et al.ll200ll ). 



In all of these PN shaping models the presence of gas heated to X-ray emitting tem- 
peratures is predicted. The advent of high-spatial-resolution, X-ray observatories (Chan- 
dra and XMM-Newton) facilitated the discovery and scrutiny of this shock-heated gas 
( jGuerrero et al.ll2005l : iKastner et al.ll2008l and references therein) and the results were sur- 
prising. The observed X-ray temperatur es of the shocked plasma have been mu ch lower, 



1 - 3 x 10 6 K, than predicted, > 10 7 K riStute fc Sahai 



2006 



i ng a variety of po tenti al solutions (s e e ISoker fc Kastnerl 120031 . for a review). lAkashi et al. 



(2006, 


2007, 


2008 


) and 


Steffen et al. 


( 


2008 



Kastner et a l. 2008 ), spawn- 



an attempt to reso lve the tempe r ature discrepancy while adequately reproducing other ob- 



served quantities. ISteffen et al.l (120081 ) investigate the role played by thermal conduction 



in est ablishing the physical properties of the X-ray emitting gas, while lAkashi et al.l (12006 



20071 ) consider the possi bility that the X-ra y emitting gas arose from a slower, post-AGB 
wind (~ 500 km/s), and lAkashi et al.l (120081 ) present numerical simulations that suggest the 
observed X-ray properties could be explained by jet-wind interaction. These newer models 
describing PN X-ray emission should, in turn, provide feedback to generalized models de- 
scribing the origin of the flows shaping PNe, which, taken together, will provide a complete 
theoretical description of the shaping and evolution of PNe. 

The identification of X-ray sources within bipolar PNe should provide unique con- 
straints on this new generation of PN shaping a nd evolution models . However, thus far 



such detections have been few and far between (IGruendl et al.l 120061 : IKastner et al.l 12001 
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Kastner el al.ll2003l ; ISahai et al.ll2003l ). During an archival program to identify serendipitous 



X-ray sources associated with PNe, we established that the XMM-Newton Serendipitious 
Surve}0 includes an X-ray source, 2XMM J174756. 2-295937, that is spatially coincident, 
within errors (~ ±2"), with Hubble 5 (Hb 5), a large, narrow- w aist bipolar planetar y nebula 
with a bright, compact, core (IPottasch fc Surendiranathll2007l ). iQuireza et al.l (120071 ) classify 
Hb 5 as Peimbert Type I PN, i.e., a He- and N-rich PN descended from a relatively massive 
progenitor star. Hb 5 is a hig h-excitation PN and, in particula r, is one of the few PNe show- 
ing the 7.652 fim Ne VI line ( IPottasch fc Surendiranathl 120071 ) . A photodissociation region 
immediat ely surrounding; the nebula may account for the presence of m olecular hydrogen 
emission ( Bernard-Salas fc Tielens 12005 ; Pottasch fc Surendiranath 2007 ). The distance to 
Hb 5 is highly uncertain. iRice et al.l (120041) derive a distance of 1. 4 kpc ± 0.3 kpc from 3D 
photoionization modeling, while IPottasch fc Surendiranathl (120071 ) argue for a much larger 
distance, ranging from 3 to 7 kpc. ICorradi &; Schwara (119931 ) find a polar expansion velocity 
of ~ 250 km s-\ while IPshmish et all (l200oh reported a faster, possibly collimated, wind of 
400 km s _1 near the center of Hb 5. 



2. Analysis 
2.1. Data 

Hb 5 lies within the field of an XMM-Newton observation of the pulsar PSR J1747-2958 
and pulsar wind nebula (MOUSE NEBULA). This observation (ObsID 0152920101) was 
performed during XMM Revolution 607 on 2003 April 02. The European Photon Imaging 
Camera (EPIC) detector arrays pn, MOS1, and MOS2 were operated in Full-Frame Mode, 
with the thick filter, for total exposures of 45.1, 51.3, and 51.3 ks, respectively. We repro- 
cessed the XMM-Newton Observation Data Files using the XMM-Newton Science Analysis 
Sofware (SAS) package version 7.1.0 with the calibration files available in Current Calibra- 
tion File Release 241 (XMM-CCF-REL-241). We filtered out high background periods and 
bad events from all observations using standard filters for imaging mode observations. The 
resulting net exposure times are 44.1, 50.6, and 50.8 ks, in the pn, MOS1, and MOS2 arrays, 
respectively. 

From the source and background regions used in the spectral extraction of the X-ray 
source at the position of Hb 5 (described in Section 2.3), we determined that the net source 
and background count rates in the energy range 0.2-2.0 keV are 1.9 ± 0.4 x 10~ 3 cnts s _1 in 



Watson ct al. 



( 20071 ) for information on the XMM-Newton Serendipitous Survey catalog. 
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pn, 8.0 ± 1.8 x 10~ 4 cnts s" 1 in MOS1, and 8.7 ± 1.9 x 10" 4 cnts s" 1 in MOS2. We merged 
the three observations using the SAS task merge. Using this merged event list, we generated 
a soft band image for the energy range 0.2-2.0 keV, binned to 5" pixels (Figure [T^,). 



2.2. Spatial Analysis 

Figure [T^, indicates that the X-ray emission at the position of Hb 5 is marginally resolved. 
Hb 5 is located ~ 7.8' off-axis in the XMM observation, and the coincident X-ray source 
displays peak emission around 1 keV. This off- axis angle and pea k emission a re wi thin the 



range of application of the XMM/EPIC PSF model described in iGhizzardil (120021 ). i.e., a 
King profile whose core radius and slope depend on energy and off-axis angle. Specifically, we 
expect a PSF core radius of 4.71" and slope of -1.44. This suggests that using a normalized 
Gaussian smoothing filter with a FWHM of 7.5" will reduce the loss of information, while 
aiding our interpretation of the detected X-ray emission. In Figure dj the smoothed X-ray 
image contours are overlaid onto the unsmoothed X-ray image (Figure [Th.), a 2MASS J-band 
image (Figure [lb), and a HST WFPC2 F658N (Ha + [Nil]) narrowband image (Figure 
[lb). The image registration was performed by the IDL Astro Librarj|^| task hastrom. A 
slight eastward shift (~ 0.45") of the HST image was required to align the stars appearing 
in both the 2MASS J-band and chip 3 of the HST WFPC2 F658N fields. Figures Hb and 
[Tb demonstrate that the peak of the X-ray emission coincides with the bright core regions 
of Hb 5 seen in the 2MASS and HST images, respectively. The apparent extension of the 
X-ray emission lies along the brightest features observed in the Ha + [Nil] image. 



2.3. Spectral Analysis 

The spatial region available for spectral extraction within the XMM image of the Hb 5 
field is severely constrained by several unavoidable artifacts. Immediately surrounding the 
region of interest are chip gaps (in all three EPIC arrays), a bright nearby soft source, and 
a readout streak from the nearby soft source. Guided by the optical HST position of Hb 5 
along with the evident soft X-ray emission, we selected an extraction region that optimized 
the signal in the spectrum while avoiding these artifacts. An annulus surrounding the source 
region and omitting nearby sources was used to estimate the background. We used the SAS 
tasks rmfgen and arfgen to generate source-specific response mat rices and e f fective area 



curves; the latter accounts for off-axis vinetting (~ 30% loss at 7.8 / . lEhle et al.l (120081 )) 



2 http:/ /idlastro. gsfc.nasa.gov 
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The resulting spectra from the three EPIC CCD detector arrays were simultaneously fit 
in XSP EC (version 12.3.1; lArnaudl ( 1996)) with a variable abundan ce thermal plasma model, 



vmekal (IMewe et al.l 



1985 



1986 



Kaastralll992l ; iLiedahl et al.lll995l ). and intervening absorp- 



tion, wa bs (IMorrison &: McCammon 



198 



We used the extinction value C 



1.60 at H/3 

2 



listed in Pottasch fc Surendiranath ( 20071 ) to fix the column density at Nh = 6.0 x I0 21 av 
and - noting the apparent presence of a strong spectral feature at ~ 1.0 keV that is likely 
due to Ne IX and Ne X lines - left as free parameters the gas t emperature and the abun - 
dance of Ne. We fixed all other abundances at their solar values !] Anders fc Grevessdll989l ). 
Although fixing the abundance of Ne to its solar value results in a global fit that is accept- 
able, this fit is poor at the spectral feature near 1.0 keV. There is a modest improvement 
in the fit when we allow the abundance of Ne to be a free parameter. This latter fit is 
presented in Figure El The temperature is found to lie in the range 2.4-3.7 MK (90% confi- 
dence range) and the Ne abundance fact or lies in the range 1 .8-6.9 (90% confidence range) 
relative to the solar abundance given by lAnders fc Grevessd (11989I ). The best model fit in- 
dicates an observed 0.2-2.0 keV flux of 7.9 x 10~ 15 ergs cm~ 2 s" 1 , and an unabsorbed flux of 
1.2 x 10~ 13 



ergs cm s . We assume a distance of 3.2 kpc after iPottasch fc Surendiranath 
J2007J) and find L x = 1.5 x 10 32 (L>/(3.2 kpc)) 2 ergs s _1 . 



3. Discussion 



Hb 5 is the fourth bi polar PNe - after N GC 7026 (IGruendl et all 120061 ) . the tran- 
sitional object NGC 7027 (IKastner et al.ll200ll ). and the possibly symbiotic system Mz 3 
( IKastner el al.l 120031 ) - with suspected diffuse X-ray emission and only the second Type I 
PN known to display X-ray emission. The other Type I PN with evidence for shocked X-ray 
emission is the [WC] PN NGC 5315, which was also disco vered serendipitously, far off-axis, 
by the Chandra X-ray Observatory (IKastner et al.l 120081 ) . Indeed, the X-ray emission re- 
gions detected in Hb 5 and NGC 5315 display similar characteristics, despite the fact that 
one is associated with a high excitation, bipolar PN and the other with a compact [WC] PN. 
Specifically, these two PNe have similar plasma temperatures (within errors) and similar 
X-ray luminosities (though the Lx values are very uncert ain, due to distance uncerta i nties) . 
In addition, both appear to show strong X-ray Ne lines. IPottasch fc Surendiranath! (120071 ) 
find an Ne abundance factor of 2.2 relative to solar in the infrared spectrum of Hb 5, similar 
to the (uncertain) overabundance tentatively determined from the X-ray spectrum. 

The bipolar PNe detected thus far have closed bipola r lobes, while many of the bipolar 
PNe observed but not detected have open bipolar lobes (IGruendl et al.l 120061 ). The role of 
open and closed lobes is an important factor that may determine key characteristics regarding 
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the presence or absence of hot gas due to wind shocks; open lobes allow the gas to expand 
and, thus, quickly cool, whereas closed lobes evidently contain t he gas, and cooling through 
heat conduction at the nebular-hot gas boundary may dominate (jSteffen et al.ll2008l ). In this 
regard, it appears significant that the extended X-ray emission in Hb 5 corresponds closely to 
the brighter regions seen in the HST image, suggesting the X-rays are coming from shocked 
material within the most tightly confined regions of the lobes. 



Recent simulations of axisymmetric jets expanding into a spherical wind by lAkashi et al. 



( 120081 ) reproduced optical features seen in some bipolar PNe. This study focused on Mz 3 and 
M 1-92, but included simulations with generic PNe parameters. Some of these key optical 
features are also seen in Hb 5, i.e., an equatorial dens e region and br i ght p olar rim. The 
evidence of a fast, possibly collimated, wind reported by lPshmish et al.l (120001 ) suggests that 
the extended X-ray emiss ion observed in Hb 5 can likely be traced to a collimated fast wind 
( lAkashi et al.ll2006l . 120071 ). The point-symmetric features of the nebula and the extension in 
the X-ray emission support the possibility of an X-ray jet or shocks from such a collimated 
flow. The spatial coincidence with the central region of Hb 5 suggests that some X-rays may 
originate from the core and/or the central star itself. However, caution is warranted, given 
the small number of counts detected due to vingetting at the large off-axis angle of Hb 5 
and the large background due to nearby bright sources convolved with the XMM PSF. 

A deep, targeted X-ray observation of Hb 5 might bring our Ne abundance factor into 
agreement with the infrared-determined value and would decrease the uncertainties in Ne 
abundance by allowing us to fit the abundances of additional elements. Also, in the archival 
XMM-Newton observation of Hb 5, small but significant (termination) portions of the bipolar 
lobes lie in the detector chip gap and the readout streak of the nearby source. Hence, a 
targeted Chandra X-ray Observatory X-ray observation of Hb 5 would provide the higher 
spatial resolution imagery that might enhance the role this PN plays in our understanding 
of the shaping and evolution of planetary nebulae. 



4. Conclusions 

Our spatial and spectral analysis of a moderately resolved, off-axis, low count XMM- 
Newton detection of X-ray emission at the location of the bipolar PN Hb 5 leads us to 
conclude that the emission originates with Hb 5. The X-ray emission is indicative of a 
thermal plasma at 2.4-3.7 MK with X-ray luminosity of 1.5 x 10 32 erg s _1 . Comparison 
with images of infrared and optical emission shows that the peak of the X-ray emission 
corresponds closely with the bright core of Hb 5 and the extension of the X-ray emission 
lies along bright extensions in the optical image. The X-ray temperature and luminosity, 
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as compared to PNe known to exhibit shock-induced X-ray emission (IKastner et al.ll2008l ) 



further suggest that the emission arises from shocks in the nebula, as opposed to the harder 



and le ss luminous X-ray emission associated with central stars (e.g. Mz 3; IKastner el al. 



( 120031 )). As only the fourth bipolar PNe with detected X-ray emission, Hb 5 may play an 
important role in the evolving study of the shaping and evolution of PNe. This PN therefore 
merits future targeted, deep, high spatial resolution X-ray observation. 
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Fig. 1. — XMM imaging of the Hb 5 field. Smoothed X-ray contours at 9, 11, 13, and 15 
counts (solid contours) and at 20, 40, and 60 counts (dotted contours) are overlaid upon 
the (a) unsmoothed XMM X-ray image, (b) 2MASS J-band image, and (c) HST F658N 
narrowband image (Ha+ [Nil]). In the lower left corner of each panel, we display the 25%, 
50%, 75% and 90% levels of the normalized Gaussian filter used to smooth the X-ray image. 
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Fig. 2. — X-ray spectra of the Hb 5 source as obtained from the three EPIC CCD detector 
arrays (crosses, diamonds, and squares; pn: top panel; MOS: bottom panel) overlaid with 
the result of the best simultaneous fit to a model consisting of a thermal plasma (vmekal) 
and intervening absorption (wabs). 



